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Abstract 

i— H Market-based instruments have been proposed as a means of implementing cost-effective and adaptive policies for biodiversity 
conservation. Their very foundation is a well-defined measure of ecological value; only the latter makes conservation services 
i 'comparable and translates them into a commodity which can be traded among market participants. But what is the right value 
when local conservation decisions affect ecological functions on a larger scale and thus create spatial externalities (site synergies) 
I on neighboring sites? By means of an agent-based model, we analyze different spatial trading rules and their implications on land 
O use decisions in a dynamic cost environment. The model contains a number of alternative submodels which differ in the individual 
5 assignment of market values and in the social organization of agents, the latter including cheap talk coordination and cooperation. 
We show that spatially explicit trading rules can effectively adapt land use towards species requirements and individual costs. For 



C/3 

a certain class of submodels (local decisions) , we find thresholds for the effectiveness of spatial incentives along a curve of critical 
. ,—i parameter values. The progression of this curve is derived analytically and shown to agree with the simulation results. Further, we 
^ show that cases of low social organization (e.g. no communication) are particularly prone to suboptimal land allocations due to 
^^misinformation and coordination failure. 

Key words: land use, spatial externalities, tradable permits, biodiversity conservation, agent based models, social interactions 
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1. Introduction 

What is the value of nature? Markets for biodiversity 
conservation are based on the possibility of determining 
standards which rate conservation services (e.g. the provi- 
sion of an acre of rainforest) in terms of their contribution 
to conservation targets. However, when ecological pro- 
cesses operate at much larger scales than typical landowner 
properties, local land use decisions are likely to produce 
side effects on the ecological value of neighboring land. This 
paper deals with the problem of including such spatial in- 
teractions into the market values of conservation measures. 

Markets are not a new choice, but still a rare one for bio- 
diversity conservation policies. Traditionally, conservation 
practice and research have been dominated by static plan- 
ning, offering a fast response to the growing environmental 
problems encompassing the progress of industrialization in 
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the last century. At present, systematic conservation plan- 
ning may include ecological interactions, external drivers 



and socio-economic conditions (Margules and Pressey 



2000 1 and is undoubtedly the cornerstone of ongoing con- 
servation efforts. A particularly important factor for the 
effectiveness of conservation planning is the inclusion of lo- 



cal conservation costs (Faith and Walker 1996 Ando et al. 



1998). Unfortunately, this insight also points to a practical 



problem of planning approaches: Local costs are frequently 
difficult or expensive to estimate, they may change over 
time, and it is seldom in the interest of landowners to 
report them honestly. In these cases, market instruments 
such as auctions, contracts or tradable permits provide an 
alternative to planning approaches because they are able 
to efficiently allocate conservation efforts, even at no infor- 
mation about local costs (information asymmetry) . 

The common denominator of market instruments is to 
introduce a measure (e.g. based on habitat quality, size) 
which rates the value of local sites or services in terms of 
their contribution to conservation targets. This measure 
translates conservation efforts into a commodity which 
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may be offered on the market ( Wu and Babcock} 1996 



Latacz-Lohmann and Van der Hamsvoort[ " 1998| Wunder 
20071. The balance of demand and supply automatically 



a) Trading procedure 



extracts local conservation costs and leads to a cost- 
effective allocation of conservation efforts as long as the 
market is competitive. 

Yet, while markets solve the problem of cost information 
asymmetry, the definition of an accurate measure of conser- 
vational value runs into problems when ecological values of 
local sites are spatially or temporally correlated, as it is fre- 



quently the case for habitat networks (Opdam et al. 2006 1 



An apparent solution is to incorporate spatial dependen- 
cies into market values. While this is generally possible, it 
implies that conservation decisions may change the market 
value of neighboring land. In the presence of such spatial 
interactions, referred to as externalities by the economic 
literature, trading might lead to inefficient outcomes if no 



correction is introduced (Mills 1980). Moreover, ecologi- 



cal processes may operate on a large range of spatial and 
temporal scales and show complex dependencies, making 



such corrections possibly very difficult. Moilanen (2005) 



for example discusses a case where static optimization is 
already computationally hard. It is therefore not clear a 
priori if market participants will be able to find their opti- 
mal allocation, in particular if they are subject to external 
drivers such as changing conservation costs. In a sense, we 
may view the process of trading conservation services as a 



self-organized reserve site selection algorithm ( |Faith et al 
2003). The aim of our analysis is to find out under which 
market rules this algorithm is successful in finding good 
solutions to the allocation problem of interacting sites, and 
under which conditions it may fail. 



2. Problem setting, tradable permits and modeling 
approach 

In this paper, we consider a tradable permit market as 
an example of a market for conservation services. Such a 
market is initiated by distributing conservation obligations 
among landowners, adding up to the desired regional level 
of conservation. Every landowner who undertakes conser- 
vation receives a certificate which documents the value of 
his contribution. Landowners who hold obligations can ei- 
ther choose to undertake conservation measures according 
to their obligation, or buy certificates from other landown- 
ers whose amount of certificates exceeds their obligation. 
Ideally, subsequent trading among market participants 
leads to a socially optimal allocation of conservation mea- 
sures across space and time while trading rules ensure that 
the conservation target is fulfilled at any time. Fig. [TJi 
shows a typical trading procedure in this kind of market. 



Translating the problem of spatial interactions to the 
context of a tradable permit market, we would like to know: 



Land use change 
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b) Spin model of local interactions 
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Fig. 1. a) Tradable permit market - Initially, sites B and E are 
conserved. E pays D to take over his conservation obligation. In the 
next step, D starts conserving his land and E can withdraw from 
conservation. This is repeated in the next step for A and B. We 
assumed here that all sites have equal value and that land use can be 
changed immediately, b) Spin model of local interactions: Conserved 
sites in direct vicinity interact with each other, leading to a higher 
ecological value than isolated sites. 



What is the right (certificate) value for a local site which 
interacts with other sites through a habitat network? Our 
approach to answer this question is based on the following 
three pillars: 

"The total amount of certificates must always reflect the 
global ecological value": Certificate trading rules must be 
chosen such that, starting from an arbitrary configuration, 
trading can only lead to a configuration of equal or higher 
value. This is ensured when the total amount of certificates 
always accurately reflects the total ecological value. 

"The individual assignment of certificates can be varied. 
The best choice is the one which realizes the desired level 
of conservation at lowest costs": Respecting the first rule, 
there might still be several ways to distribute the total 
amount of certificates among local landowners. The indi- 
vidual assignment should be chosen to maximize the cost- 
effectiveness of the land allocation emerging from trading. 

"Different social structures might affect the optimal choice 
of assignment rules": As landowners' decisions are mutu- 
ally dependent, aspects such as communication may play 
an important role for individual decision making. Different 
social organization structures may suggest different indi- 
vidual assignment rules, even if the ecological situation is 
the same. Fig. [2] illustrates how we translated these ques- 
tions into our simulation. In the rest of this section, we will 
give a non-technical motivation for the chosen submodels. 
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Fig. 2. Modeling approach: On the left the ecological evaluation, 
starting from a land configuration with local sites A-E and leading 
to an estimate of ecological value of the whole configuration. On the 
right the market, starting from certificate values for the sites A-E and 
leading to a changed land configuration. Both sides are connected 
by the individual assignment. 



2.1. Submodel 1: Ecological evaluation 

The survival of species in fragmented, dynamic land- 
scapes is subject to a complicated interplay of reproduction, 
dispersal and mutual interaction. To estimate the value of 
such landscapes, we need to define categories in which we 
rate them. These categories are by their nature norma- 
tive and subject to the current understanding of ecosystem 
functions and societal preferences. An extensive amount of 



literature deals with this question, e.g. ( 


Saunders et al. 


1991 ||Frank and Wissell|1998||Opdam et al. 2003 Brugge- 


man et al. 


2005 


). Here, it is widely agreed that habitat 



connectivity is a key factor for the viability of a large range 
of species. For our simulation, we will make a crude sim- 
plification: We assume connectivity and area are the only 
important factors for the ecological value of a site. More 
explicitly, we assume that N sites of unit size, assembled to 
a network, have an ecological value of 
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2.2. Submodel 2: Individual assignment 

Having defined a value U for any landscape configura- 
tion, we now ask: How should U be split up among local 
sites? We consider two submodels which distribute the 
same total amount of certificates, but differ in the individ- 
ual assignment of certificates among agents: 

Additive certificate values (1): A straightforward way is 
distributing any gain in U that results from mutual bene- 
fits of neighboring sites among the neighbors involved (See 
Fig. [3]) . The proportion assigned to each landowner should 
reflect the contribution of his land to the ecological value 
of neighboring sites. It is easier to assume that mutual ben- 
efits are symmetric between sites, but the non-symmetric 



case (see Vuilleumier and Possingham 2006) can equally 



be treated. In the context of payment schemes, a similar 



rule has been proposed by Parkhurst et al. ( 2002 ) . We call 
this option " additive certificates" because the value of all 
distributed certificates adds up to the global value and 
the value of any subset of sites is the sum of its certificate 
values. 
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(2) 



Marginal certificate values (2): For the second option, we 
set the certificate value of a site i equal to the marginal 
change of global value resulting from a change at this site 
while all other sites remain as they were: 



U' 



old i 



(3) 



Unlike additive certificates, the marginal rule assigns all 
costs and benefits created by a land use change to its orig- 
inator. Yet, while such a situation is generally favored by 
economic theory, practically this is connected with consid- 
erable complications: Unlike additive certificates, marginal 
values are in general not additive. This means that the 
marginal value of 2 sites is not generally the sum of the 
two sites' marginal values. Further, the resulting payoffs for 
landowners depend on the order of trading. This does not 
pose theoretical problems, but it might lead to severe ac- 
ceptance problems among landowners. Both properties are 
illustrated in Fig. [3] 



(1) 2.3. Submodel 3: Social Organization 



where Pi is a parameter measuring the connectivity of the 
i-th site, while to, taking values from to 1, measures how 
much the species under consideration values habitat con- 
nectivity over area. Despite its simplicity, this function re- 
flects the importance of area and connectivity as an essen- 
tial insight from spatial ecology. Qualitatively, our results 
can be expected also to hold true for more complex evalua- 
tion functions such as the one derived by Frank and Wissel 
(20021. 



Clearly, if certificate values of a site depend on the neigh- 
boring sites, it is important for landowners to know their 
neighbors' intentions. But how much does social organiza- 
tion affect the results of the simulation? We assume that 
agents are always profit maximizing and myopic in the 
sense that they base their decision on the best action for 
the next timestep without displaying strategic behavior. 
Within this setting, we consider three behavioral submod- 
els: 
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Fig. 3. Two sites A and B of equal size pose a mutual benefit to each 
other. When isolated, both sites have an ecological value of 2 each. 
Connecting the two sites creates an additional value of 2, resulting 
from a benefit of 1 in each direction. Given additive certificate values, 
each landowner receives 3. The marginal value of a conserved site in 
the absence of any other site is only 2, while the second site adds all 
mutual benefits and thus receives 4 — the order of creation matters. 

Null model (a): In the null model, agents observe the 
present landscape configuration and decide in the prospect 
of the future land configuration being the same as the 
present one. 

Coordination (b): It is natural to assume that agents will 
communicate their future intentions. This is beneficial 
because it increases the accuracy of the estimate about 
the upcoming land configuration. Coordination is under- 
stood by us as communication of non-binding information 
about the present state of the decision, often called cheap 



talk (Farrcll and Rabin 19961. Experimental studies have 



shown that the possibility of coordination by cheap talk 
leads to an increased probability of finding cost-effective 
configurations (Parkhurst et al. 20021. In general, cheap 



talk also offers the possibility of strategic lies. This option 
is omitted in the simulation and hence agents will always 
stick to the action they communicated as long as condi- 
tions do not change. 



Table [T] summarizes the discussed model options and their 
properties. For cooperation, there is no difference between 
the two certification processes, because these differ only 
in the individual assignment of certificates, but not in the 
total assignment. 

3. Model description 

3.1. Overview and purpose 

To represent interactions between site values, we use a 
2-dimensional spin model with local interactions (see Fig. 
[lj>). Spin models have been used repeatedly in the litera- 



ture to analyze phenomena of social interactions (Galam 



and Zucker 



ct al 



2000 



[2000] |Sznajd-Weron and Sznajd| [20001 



Holyst 



Schweitzer et al. 1 12002 ) . The model we present 



differs from a standard Ising model in two major points: 
Firstly, its dynamics are driven by a randomly fluctuat- 
ing external field at each local site, representing changing 
conservation costs. Secondly, the system is constrained by 
the conservation target. To our knowledge, there has been 
no systematic analysis of this model in the literature, al- 
though it exhibits similarities to the Random Field Ising 



Model (Imry and Ma 19751 and non-equilibrium models 



such as Hausmann and Rujan ( 1997 ) and |Acharyya ( 1998 1. 
However, our focus is not to solve a somewhat modified 
Ising model. Instead, technically we are comparing differ- 
ent algorithms (specified by certification rules and agents' 
behavior) in terms of costs they require to reach the con- 
servation target. 



Symbol Connotation 



Range 



i-th cell on the grid 
State of the i-th cell {0, 1} 
Costs of (Tj = 1 at t [1 - 8...1 + <5] 
Market price [1 - 6...1 + S] 



Ci(t) 

P 



Cooperation (c) : Further payoff improvements are possible 
if agents did not only coordinate, but cooperate. Coopera- 
tion means that conservation is provided if it is beneficial 
for the group, even when it is not in the immediate interest 
of the individual. In such a system, agents would reveal 
their true costs to the group, the best configuration would 
be chosen, and the payoffs would be distributed among 
the group members. However, this creates the temptation 
to communicate higher than the true costs to increase the 
individual share of the group benefits. Cooperation is a 
public good game and cooperators have to fear exploita- 
tion by defectors and free-riders. We will address this issue 
in detail in the discussion. 

For the first two options a) and b), decisions are made 
locally, while cooperating agents c) decide on a global level. 
We will see later that these structurally different decision 
processes also reflects in the resulting land use pattern. 



8 Cost heterogeneity [0..1] 

m Connectivity weight [0..1] 

A Ecological target [0..f] 

Table 2 

List of state variables (top) and parameters (bottom) 

3.2. State variables and scales 

The simulation is carried out on a 2-dimensional grid 
with 50x50 grid cells and periodic boundary conditions. 
We also refer to the N=2500 grid cells as sites. Every grid 
cell Xi is owned by a different landowner (agent) and can 
be occupied with a habitat (cr, = 1) or be used for other 
purposes (<Xj = 0). The occupancy of a grid cell results 
in conservation costs Ci(t) which have to be borne by its 
landowner. These costs are sampled at each time step from 
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Label Submodel 



Decisions Game Type Path Communication 



(a.l) Null additive 

(a. 2) Null marginal 

(b.l) Coordination additive 

(b.2) Coordination marginal 



1-player 



independent 



None 



Local 



dependent 
independen 

i 

dependent 



independent 
Coordination Intentions 



(c) Cooperation 



Global Public Good independent Costs 



Table 1 

Overview of alternative submodels for the certification processes and agent behavior 



a uniform distribution of mean 1 and width 5. Following eq. 
[T] we assume that the conservational value of a landscape 
configuration, given by the values of (<7i . . . cr/v), is 



N 

U = (1 - m) °i 



N 

E 

i,<3>, 



(4) 



where < j >i indicates all j which belong to the 8 cells Xj in 
a Moore neighborhood of Xi and m weights the importance 
of connectedness over area for the relevant species. The 
market is constrained by a conservation target 



U = X-N . 



(5) 



which is set by the parameter A. Landowners demand cer- 
tificates until they reach the target A • N and thereby cre- 
ate an equilibrium price P per unit value. A list of all ba- 
sic state variables and parameters is given in Table [2] For 
convenience, we define the following abbreviations for the 
average number of occupied sites and the occupancy of the 
neighborhood of Xi, respectively: 



1 N 1 

i=l 



(6) 



<3>i 



Further, we define connectivity K as the average neighbor- 
hood occupancy (3 of the occupied sites 



K 



1 N 



a-N 



(7) 



i=i 



and societal costs C of a landscape configuration as the 
sum of the costs of all conserved grid cells, divided by the 
number of grid cells. 



1 N 



(8) 



3.3. Process overview and scheduling: 

At each time step, new costs Ci are sampled. Agents de- 
cide to maintain a site as habitat based on their costs Cj, 
the market price P and the estimated certificate value bi. 
They maintain a habitat on Xi at timestep t if 



7r» ee -a(t) +P -bi(t) > 0. 



As discussed in section |2.2[ we consider two certification 
options. Assuming symmetric benefits, eq. [4] results in ad- 
ditive certificate values of 



bi = (1 — to) + m ■ Pi . 



(10) 



The marginal certificate value of a site can be calculated 
from eq. [4] according to eq. [3] 



b[ = (1 - to) + 2 • ft 



(11) 



Although the process of certificate exchange is not mod- 
elled explicitly, the reader may assure himself that eqs. 
10 and [iT] do exactly create the decision base in terms of 
costs and benefits to match the trading rules described 
in section |2.2| in particular also the path-dependence of 
marginal certificates. The final benefits bi rewarded at the 
end of the round can differ from the estimated benefits bi 
because subsequent decisions by other agents can change 
the landscape configuration. The accuracy of the estimate 
bi depends on agents' communication and organization 
abilities. The three submodels discussed in section 1231 are 
implemented in the following way: 

Null model: All agents decide in parallel according to eq.|9] 
without being informed about the decisions of other agents 
at this timestep. 

Coordination: Agents decide sequentially in random or- 
der according to eq. [9] After each decision, all agents are 
informed about the new configuration. This procedure is 
repeated a number of times, mimicking the outcome of a 
non-binding exchange of information on the basis of will- 
ingness. 

Cooperation: Cooperation is modeled by global optimiza- 
tion, which is described in the next subsection. It is as- 
sumed that all agents communicate their costs and decide 
together. 

Step by step, the decisions of all agents are collected and 
the resulting ecological value U as given in eq. [i] is calcu- 
lated. U is constrained by the conservation target of A • N 
(eq. [5]) . The emergence of an equilibrium between fixed de- 
mand (target constraint) and supply (eq. [9]) is modeled by 
adjusting the market price P until the agents as a whole 
choose to conserve the target value A • N (Drechsler and 



Watzold 20071. The order in which agents are asked is 



(9) randomized at every timestep, but does not change while 



5 



the market price is adjusted. Fig. [4] shows a flow diagram 
of the processes within one timestep. 
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Fig. 4. Flow diagram of processes within one timestep. 

3.4. Initialization and Stochasticity 

Simulation runs were initialized with a random landscape 
configuration on the target level. A series of tests showed 
that the initial configuration has no influence whatsoever 
on the resulting landscape after a sufficient amount of time. 
This holds also true for non-random start configurations as 
well as random initial ccopoint levels. Data acquisition was 
started after 600 trading steps to ensure that the simulation 
had reached its steady state. The results show the mean of 
50 runs. Standard deviations for all values were calculated 
but omitted for presentation because they were small and 
showed no significant information. 

3.5. Optimization 

To model a perfectly rational cooperating group of 
agents, we performed a global optimization of the land 
configuration. We used a slightly modified simulated an- 
nealing algorithm, which delivered better results than the 



original algorithm by |Kirkpatrick et al. (19831. The goal 



of the optimization was to minimize the societal costs 
as defined in eq. [8] under the constraint of meeting the 
regulatory target as given in eq. [4] 




U=X-N 



A random site Xi was occupied with probability 
p{(Ti,T) = mm 1 1, e T 



(12) 



(13) 



ration. The simulated annealing was performed in n steps 
with a new random subset of N ■ v of the sites at each step. 
To satisfy the constraints, each step was followed by adding 
(removing) conserved sites starting with the sites of high- 
est (lowest) c/> until the target value of U — A • N was met. 
Temperature decay was exponential with decay parameter 
r per time step. Table [3] lists the parameter values used. 



Parameter Connotation 



Value 



where 4>i(ai) is the marginal benefit cost ratio of site Xi and 
$ is the average benefit cost ratio of the present configu- 



Trj Start temperature 3 

r Decay parameter 0.001 

n Steps 3500 

v Fraction of sites 0.5 

Table 3 

Optimization parameters and standard values chosen for the simu- 
lated annealing 

4. Analytical Bounds Clusters and Disorder 

Although the number of landscape configurations which 
fulfill the conservation target is extremely large, two struc- 
tures occupy an outstanding position: One is a landscape 
where habitats are concentrated into one big cluster, the 
other where habitats are scattered randomly, disordered 
according to the lowest costs. Examples of these states will 
appear later in the results (Fig.j5|. Clustering and disorder 
mark the two outer boundaries in terms of possible connec- 
tivity values: No other state produces more connectivity 
than a cluster, and no reasonable state (taking aside anti- 
correlated structures) produces less connectivity than the 
scattered, disordered state. It will prove useful for the inter- 
pretation of the results to derive some analytical properties 
of these two states. 

4.1. Critical Values 

Let us assume that agents can only choose between clus- 
tering and spread. Clustered structures lead to a higher 
certificate value per cell, but also to higher average costs 
per cell, because a spread, disordered configuration can 
more effectively concentrate conservation on the sites with 
the lowest costs. At low values of cost heterogeneity r5 com- 
pared to connectivity weight m, a clustered structure is 
clearly favored. At increasing cost heterogeneity, we expect 
a critical value S c where the net benefits from clustering 
become smaller than the net benefits from spread. 

We can derive this critical value by equating the worst 
benefit-cost ratio of cells within a cluster with the best 
ratio of an isolated cell. In a cluster, the habitats with the 
highest costs have c = 1 + 6, while outside the cluster, the 
cells with lowest cost have c = 1 — 5. The ecological value of 



a clustered and a disordered cell is given by eqs. 10 and 11 
Hence, we have 



G 



1 

1 + m 
1 + *S 



1 — m 
1 — r$ c 

1 "add 
1 — 771 
l-6 c 



<^add = m/(2 - 777) 



(14) 



(15) 



as critical values for additive and marginal certificates, re- 
spectively. 

4.2. Clustered and disordered cost level 

Further, we are be interested in the costs of maintaining 
the land at either of the two states. Let c p (Q denote the 
inverse cumulative distribution of the chosen probability 
density function of local costs. In economic terms, this can 
be regarded as the marginal cost function c p (£) of conser- 
vation. With costs being distributed uniformly across the 
interval [1 — S ... 1 + 5), c p (£) is given by 



cp(<) = (1 - S) + 26( . 



Cells in a cluster have 8 neighbors and therefore yield an 
ecological value of 1 per cell (eq.[4}. Therefore, a number of 
A • N patches meets the regulatory target of an ecological 
value of A per cell (eq. [5]). As costs are uncorrelated, the 
mean costs c p within a cluster are approximately equal to 
the mean costs of the landscape (for the chosen function 
c p — 1) as long as finite size effects can be neglected. Thus, 
the societal costs of reaching the target of a value A per cell 
by means of a cluster are 



Cc 



1 

N 



c p (()d(-N\ = A 



(17) 



In the disordered state of density a, where occupied sites 
are distributed randomly according to the lowest costs, each 
occupied cell has on average 8 • a neighbors, leading to an 
average certificate value of (1 — to) + to • a per grid cell. To 
reach a target of A per cell, we require that 



a {(1 — to) + to ■ a} = A 
which leads to 



-1 + m + J (-1 + to) 2 + 4toA 

a = . 

2to 

The societal costs in the disordered state become 



(18) 



(19) 



C do = aJ c p (()d( = a-{l-(l-a)-S) (20) 
o 

with a given in eq.[l9j The two cost curves for the clustered 
and the disordered state mark an upper boundary for cost- 
effective configurations. As a function of S, they intersect 
at 

5=^—4= 777 (21) 

which coincides with the critical point for marginal certifi- 



5. Simulation Results 



5.1. Critical Values 



Scanning the parameter space of to and 6, we observe 
steep transitions of all aggregated state variables along a 
curve of critical values S c when moving from large to low 
ratios of 777, over <5. Parameter values beyond this curve 
lead to disordered landscape structures, while values be- 
low 5 C lead to ordered, connected structures. As expected, 
the shape of the transition line differs for additive and 
marginal certificate values. Fig. [5] shows the simulation 
results for additive (a.l) and marginal (a. 2) certificates 
at A = 0.1 together with the analytically derived curves 
eqs. [14] and [15J The right side of Fig. [5] shows a vertical 
cross-section of the latter plots in 5-direction at to = 0.5. 
The transition values 6 C in this curve agree with the theo- 

and [15] (Table @|. 



(16) retical values from eqs. 



14 



The introduction of coordination (submodel b.l and b.2) 
leaves these results largely unchanged. Fig. [6] shows that 
coordination only slightly increases the tail after the tran- 
sition point, but does not change the transition point. In 
contrast to the steep transitions of local decision processes, 
submodel (c) with cooperative decisions shows a rather 
linear transition for the same parameter values. 



Model 



Simulation Theory 



Additive Value 0.33 ± 0.02 0.33 
Marginal Value 0.50 ± 0.02 0.5 

Table 4 

Critical value S c (measured as <5 at half transition) at m = 0.5 from 
Fig [5] together with theoretical expectations 
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Fig. 6. Connectivity (on the y-axis) as a function of cost heterogeneity 
8 at a fixed m = 0.5. It can be seen that increased communication 
does not affect the critical point <5 C of the additive certificate curves. 



cate values (eq. 15) 
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Fig. 5. On the left, landscape connectivity for marginal certificates (top), additive certificates (middle) and cooperation (bottom) as a function 
of cost heterogeneity 5 and connectivity weight m at A = 0.1. Darker colors indicate lower connectivity. Dotted lines display the theoretical 
curves as given in eq. 1 1 4| The right graph shows a cross-section of the three plots in 8 direction at m = 0.5 together with typical landscape 
structures emerging from the simulations in the three domains: Clusters (left), transition states (middle) and disorder (right). 

5.2. Cost-effectiveness 



One major question was how the three options would 
perform in terms of costs necessary to reach the conser- 
vation target. Figure [7] shows the societal costs as a func- 
tion of the cost heterogeneity a at m = 0.5, together with 
the theoretical cost levels for the clustered and the disor- 



dered state as calculated in eqs. 17 and |20| (dotted lines). 
Noticeably, the cost function of additive certificates dis- 
plays a hump around the transition area (Fig. [7]), leading to 
approximately 20% higher costs for additive compared to 
marginal certificates. When considering the theoretical cost 
functions for the ordered and the disordered state (eqs. 17 
20 ), the reason becomes evident: Marginal certificate values 
switch from clustering to disorder right at the intersection 
of the cost lines (eq. 21 1. In contrast, additive certificates 
switch already at a lower cost heterogeneity S. Here, the 
cost level of a disordered configuration is still higher than 
of a clustered one, resulting in efficiency losses in this area. 
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Fig. 7. Societal costs as a function of cost heterogeneity <5 at m = 0.5 
and A = 0.1 (solid lines) together with theoretical costs for the 
clustered and the disordered state (dashed lines). 



5.3. Coordination and cooperation 

The inefficiency which has been observed for additive 
certificates is to a great extent mitigated by the introduc- 
tion of coordination (Fig [8]). For a value of 10 commu- 
nication steps per trading period, only small differences 



between additive and marginal certificates remain. This 
results mainly from a better adaptation of land use to the 
current costs. Transition points and thus the landscape 
structure remain largely unchanged by coordination as 
discussed on Fig. [6] 

Not surprisingly, the societal costs under cooperation are 
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considerably lower than for the other options (Fig. j7|. 
They differ mostly around the critical values S c where local 
decisions (submodel a and b) are more difficult than for 
extreme values of S where one allocation pattern (cluster 
or disorder) is clearly favored. 



0.13 
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Fig. 8. Effect of coordination: The societal costs are plotted as a 
function of cost heterogeneity 8 at m = 0.5 and A = 0.1. With 
increasing coordination steps, additive certificate costs approach the 
costs of marginal certificates. 



simulation outcomes merely to the clustered or the disor- 
dered state as described above. In contrast, simulations 
of coordinating agents (submodel c), as an example of a 
global decision process, generate the whole range of in- 
termediate clustering throughout the parameter space. 
Although not presented here, we checked that temporal 
and spatio-temporal correlation of costs generated by a 
random walk algorithm may also smooth out the steepness 
of the transitions in favor of more linear shapes. The rea- 
son is that when cost correlation length is large compared 
to the ecological correlation length, the chance of getting 
stuck in a local minimum is rapidly decreasing. 

Cost differences between the different submodels arise 
most evidently in the transition regime, where connected- 
ness becomes just as beneficial as spread. Here, additive 
certificates lead to a suboptimal level of landscape connec- 
tivity and result in increased societal costs. The origin of 
theses cost differences lies in the different spatial incen- 
tives of both options: Additive certificates assign a part of 
the spatial benefits of a site to its neighbors and therefore 
give less incentive for clustering from the individual view- 
point. Yet, the strong costs increase can be predominately 
explained by losses due to misinformation. At increased 
rates of coordination, cost differences between additive 
and marginal certificates decrease strongly. 



6. Discussion 

6.1. Main findings 

We have introduced a simple model of a market for 
tradable conservation services with spatially correlated 
certificate values. The model enabled us to compare differ- 
ent spatial certification rules such as additive and marginal 
certificate values and different types of agent behavior, in- 
cluding local and global decision processes. We found that 
local decision structures (null model a.l,a.2 and coordina- 
tion b.l,b.2) lead to a partition of the parameter space in 
two areas with very distinct landscape patterns, divided 
by a small transition area. These predominant states are 
a completely clustered state, appearing when cost hetero- 
geneity S is low compared to the connectivity weight m, 
and a completely disordered state, emerging when <5 is high 
compared to m. 

The reason for the stability of these states in a dynamical 
and heterogeneous cost background lies in the fact that, 
starting from one of them, agents repeatedly end up in 
a very similar configuration since trading stops as soon 
as a local cost optimum is found. Even in the case of co- 
ordination with unlimited communication steps, agents 
stop trading once trade cannot improve the position of 
any single agent (Fig. |6|. This limits the possibility of 
finding more complicated allocation patterns in a highly 
stochastic environment and explains the confinement of 



A more efficient allocation of conserved sites can only 
be reached through cooperation where the configuration 
space is explored in a way which avoids getting stuck in 
local minima. However, this requires landowners to reveal 
their true costs and temporarily accept lower payoffs if the 
group payoff is increased in exchange. Payments between 
landowners (side payments) could compensate for these 
losses, but it is an open question how this kind of coopera- 
tive system should be organized and if it is stable against 
exploitation by defectors and free-riders. Moreover, it is 
not clear whether defection can be detected. If agents can 
estimate the cost of neighboring landowners, defection is 
directly detectable and it seems reasonable that mech- 
anisms such as reputation and social norms are strong 



enough to facilitate cooperation ( 


Fehr and Schmidt 


1999 


Sigmund et al. 


2001 


Milinski et al. 


2002). However, 


if de- 



fection is only indirectly detectable, much weaker support 



is present for cooperation to hold (Sell and Wilson 1991 1. 



Concluding, we believe that the general problem of in- 
cluding spatial interactions into a local value function is 
largely independent of the exact shape of the ecological 
benefit function. Consequently, these results will qualita- 
tively hold true for a wide range of ecological systems with 
spatial interactions. However, one dimension has been ne- 
glected in this simulation: Time, both on the ecological 
side (turnover, recreation times of habitat) as well as on 
the economic side (time preferences of humans) will be 
of major importance for a structurally realistic model of 
a conservation market. An extension of the model to a 
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market which includes space and time for evaluating con- 
servational value should therefore be considered for future 
research. 

6.2. Policy implications 

Tradable permit markets for biodiversity present a 
promising policy instrument when local conservation cost 
are not known and flexibility and adaptation are of im- 
portance. However, it is too short-sighted to assume that 
giving the right exchange rate to a service such as bio- 
diversity conservation will automatically lead to efficient 
allocations of conservation efforts. To find robust market 
designs, what is needed are realistic simulations incorpo- 
rating both ecological and socio-economic aspects down 
to the individual landowner level at which important pro- 
cesses take place. Although we have shown that marginal 
certificates perform better than additive ones, considera- 
tions which are beyond the scope of this paper might still 
favor the latter in practice. As an example, marginal ben- 
efits might be considered unfair by landowners because 
payoffs depend on the order of trading. In the end, markets 
must not be judged by the standards of a first-best world, 
but by practice in the field, where other policy options also 
cause severe efficiency problems. If beneficial social struc- 
tures are present, such as a strong tradition of cooperation 
and good management capacities, markets for biodiversity 
may be superior to other instruments, while under differ- 
ent circumstances, other policies like taxes may be more 
successful. 
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